Charge Transport Processes in a Superconducting Single-Electron Transistor Coupled 

to a Microstrip Transmission Line 
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We have investigated charge transport processes in a superconducting single-electron transistor 
(S-SET) fabricated in close proximity to a two-dimensional electron gas (2DEG) in a GaAs/AlGaAs 
heterostructure. The macroscopic bonding pads of the S-SET along with the 2DEG form a microstrip 
transmission line. We observe a variety of current-carrying cycles in the S-SET which we attribute 
to simultaneous tunneling of Cooper pairs and emission of photons into the microstrip. We find 
good agreement between these experimental results and simulations including both photon emission 
and photon-assisted tunneling due to the electromagnetic environment. 
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Charge transport in all-superconducting single- 
electron transistors (S-SETs) has been studied in 
detail, and numerous current-carrying cycles have been 
identiE|ed.|^including the Jusephson-quasiparticle (JQP|)| 
cycle,EJl3'LlQ the 3e cycleB and Andreev reflection.El 
What processes occur in a given S-SET is determined 
by a complex interplay between Cooper pair tunnel- 
ing and charging effects, leading to rich behavior in 
S-SET current-voltage characteristics. Further com- 
plexity is added by the electromagnetic environment, 
which can cajise photon-assisted tunneling (PAT) of 
quasigiarticlegim and incoherent tunneling of Cooper 
pairs.Bll2l Experiments have addressed the influence of 
the environment on a varietit-pf tunnel junction systenis, 
including electron pumpsptil single tunnel junctions,t2l 
Josephson junction arraysll^ and S-SETs fabricated- 
above a ground plane, focusing on the supercurrent.til 
The effect of the environment on S-SETs remains 
interesting given the potential use of-Coulomb blockade 
devices for quantum computation.t£l Here we report 
measurements on samples similar to that in Ref. 
but focus on transport at higher biases. We observe 
previously unreported charge transport cycles which we 
attribute to Cooper pair tunneling with simultaneous 
emission of a photon into the environment. We find good 
agreement between our measurements and simulatiopa 
based on the orthodox theory of electron tunneling Jlj'tZl 
including the effects of both photon emission and PAT. 

Our samples consist of Al/A10a;/Al S-SETs fabricated 
on a GaAs/AlGaAs heterostructure containing a two- 
dimensional electron gas (2DEG) located ~ 50 nm below 
the substrate surface, as illustrated in Fig. |l|. The het- 
erostructure, grown on a GaAs substrate using molecular 
beam epitaxy, consists of the following layers: 1000 nm 
of GaAs, 47 nm of Alo.3Gao.7As and 5 nm of GaAs. The 
Alo.3Gao.7As is delta-doped with Si 22 nm from the lower 
GaAs/Alo.3Gao.7As interface, at which the 2DEG forms. 
The SETs are fabricated using electron-beam lithography 
and shadow evaporation techniques. The central island 
of the SET is coupled by a capacitance C2DEG to the 
2DEG beneath it, which is held at ground. The SIi3j is 
surrounded by six Au gates, as described elsewhere.EJ 



We model the SET as shown, with tunnel junction re- 
sistances and capacitances i?i(2) and ^1(2) respectively. 
A voltage Vg applied to one of the Au gates is used to ad- 
just the electrostatic potential of the SET island through 
the capacitance Cg. The remaining gates are held at 
ground and contribute Cs to the total island capacitance 
= Ci + C2 + Cg + Cs + C2DEG- We focus first on the 
2DEG far from the SET, which in conjunction with the 
macroscopic bonding pads forms a microstrip transmis- 
sion line to which the SET is coupled. Such microstrip 
can resonate at microwave frequencies .and .thereby affect 
transport in tunnel junction systems.E3li3 We have ex- 
amined two such samples with junction resistances rang- 
ing from 24 kfl to 192 kfi, and a third control sample 
fabricated on a piece of GaAs without a 2DEG. Sample 
parameters are given in Table ^ below. 

The large value of Czdeg « 356 aF (382 aF) for 
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FIG. 1: Schematic diagram of a sample, showing the macro- 
scopic bonding pads and Au gates. Expanded view: cross- 
section of the tunnel junctions and the capacitance C2DEG 
coupling the SET to the 2DEG. Circuit diagram: specific 
parameters discussed in the text, including the parallel com- 
bination of Cs and C2dbg- 
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FIG. 2: (a) Measured / for Sample 1 on a logarithmic scale 
vs. V and Qo- (b) Simulated / using the parameters in Ta- 
ble |, without PAT. (c) Simulated / including PAT, using 
-Re = 9 n and Te = 3 K. Lines correspond to quasiparti- 
cle (solid), Cooper pair (dotted). Cooper pair with photon 
emission (dashed) and 3e (dash-dotted) transitions. 



Sample 1 (2), along with the total gate capacitance 
Cs + Cg ^ 20 aF helped ensure that the charging energy 
Ec — e-^/2Cs satisfied Ec < A where A is the super- 
conducting gap. For the control sample, Ci and C2 were 
sufficiently large to give a comparable Ec- For all the 
samples, the Josephson energy given by the Ambegaokar- 
Baratoff relationtj satisfied Ej. = hA/8e^Ri < Ec so 
that the charge state n of the island is well defined, jus- 
tifying use of the orthodox theory. We define the charg- 
ing energy of the SET U{n) = Ec{n — Qo/e)^, where 



TABLE I: Physical parameters of the samples. Energies are 
in i-LeV, capacitances in aF and resistances in kQ. 



Sample 


A 


Ci 


C2 




-Ri + R2 


R1/R2 


1 


207 


181 


120 


118 


305 


0.59 


2 


201 


375 


260 


77 


63.5 


0.6 


control 


198 


515 


322 


91 


91. 


0.7 



Qo = CgVg is the offset charge induced by the gate volt- 
age. At low temperatures, electrons will tunnel through 
a given junction only when the change in free energy 
Sf — ff — fi > 0, including both changes in U(n) and 
work performed by the bias voltage V. For V < 4A 
single-electron tunneling is in general forbidden by the 
combined effects of the Coulomb blockade and the super- 
conducting gap, though more complex charge transport 
processes such as the JQP and 3e cycles are allowed at 
specific locations in the Qq-V plane. 

We measured the samples in a dilution refrigerator at 
Tlnix — 20 mK using a four probe configuration symmet- 
rically biased with respect to ground. High-frequency 
noise was excluded, using standard techniques as de- 
scribed elsewhere.EJ Sample parameters such as the ca- 
pacitances and summed junction resistances i?i +R2 were 
determined from the measurements. The resistance ra- 
tio R1/R2 was chosen to give the correct relative size of 
features associated with each junction in the simulations. 
Current / was measured versus bias voltage V for a se- 
ries of different values of the offset charge Qo for all three 
samples; results for Sample 1 are presented in Fig. ||(a). 

Several resonances are visible for voltages between 500 
and 800 fiV. The most pronounced are the usual JQP 
peaks, some of which are paralleled by two smaller res- 
onances. We believe these additional resonances are 
analogs of the usual JQP peaks, but differ in that Cooper 
pairs emit one or more photons into the transmission line 
while tunneling. We refer to them as JQP-ph peaks. Sev- 
eral smaller features are also visible near Qo — 0-5 at 
biases between 300 and 450 /iV. To support our inter- 
pretation of the JQP-ph resonances, and to identify the 
smaller features, we begin with conservation of energy, in- 
cluding not only the usual terms accounting for cjjeation 
of quasiparticles and changes in charging energyjjcj but 
also a term allowing for emission of photons: 

KiTOicy = U{n + Sm) — U{n) + qA + khujs (1) 
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Here Kt = 1 — (Ci + Cg)/2 is the fraction of the bias 
voltage available across junction i and nii the number 
of electrons transferred across it, dm = 7712 — toi is the 
change in the island charge, q is the number of quasipar- 
ticles created, and k is the number of photons of energy 
hujs emitted. The proper choice of mi, TO2, q and k, gives 
transition rules for quasiparticle. Cooper pair. Cooper 
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FIG. 3: Tunneling processes discussed in the text, for posi- 
tive bias voltage and increasing island charge. 



pair with photon emission, and 3e tunneling processes. 
These are indicated schematically in Fig. ^, and transi- 
tion lines corresponding to them are shown in Fig. ^. 

The JQP features in Fig. ||(a) lie along Cooper pair 
transition lines, either with (J-ph) or without (J) pho- 
ton emission, since Cooper pair tunneling is a resonant 
process. For the primary {k — 0) JQP peak between 
Qo « —1.0 and —0.6, the cycle begins with a Cooper 
pair transition through Jl (—1 —* —3), followed by suc- 
cessive quasiparticle transitions through J2 (—3 — > —2 
and —2 —1) which return the island to its original 
charge state and allow the cycle to begin again. As it 
must, the JQP peak begins just above the quasiparticle 
line for the —2 —1 transition, which gives the require- 
ment eV > {Ec + 2A) w 533 /ieV. 

To address the JQP-ph peaks, the fundamental fre- 
quency for photon emission was determined by adjusting 
hujs so that a one-photon J-ph line coincided with the 
first JQP-ph resonance between Qo = —1.2 and —0.8. 
Best results were obtained for huJs = 136 fieV. For 
Cooper pair tunneling with fc = 1 we obtain the condition 
eV > [Ec + 2A + fiLOs/'^.) ~ 601 fieV, in good agreement 
with the observed location of the resonance. Since the 
typical lead width (> 20 /im) for the pad structures is 
much greater than the 2DEG depth of 50 nm, the-jcf- 
fective wavelength for excitations on the microstripEH is 
Ae ~ 2Trc/^/eoJs « 2.5 mm where e ~ 13 is the dielec- 
tric constant of GaAs. The fundamental resonance cor- 
responds to a length £ = Ae/2 « 1.25 mm, in agreement 
with the known dimensions of our lead structures. Reso- 
nances with a best fit for Tiujs = 160 /ieV were observed 
in Sample 2, while none were observed in the control 
sample, supporting our interpretation of the data. The 
position of the k = 2 J-ph line does not match well to 
that of the associated resonance, perhaps not surprising 
given its small size and the large background current. 

To gain further insight into the features in Fig. |l](a), we 
turn to our simulations. Here, to account for the JQP-ph 
resonances, we make a phenomenolagical modification of 
the theory of Averin and AleshkinE3n3 and write for the 



cp 



2 



^q.pEj 



kfiuJsY 



(2) 



where 5f is the change in free energy due to the tunnel- 
ing of the Cooper pair and Fqp is the rate of the first 
subsequent quasiparticle tunneling. The Ak were cho- 
sen as ^0 = 1 (giving the usual Averin-Aleshkin lerm) 
and Ak = ^{Ec/fiujs)'' eyi]i{—Ec/'hu}s)^ the weightslliior 
independent emission of photons into a single mode.EJ 

Results of a simulation including the effects of pho- 
ton emission are shown in Fig. ^(b). The agreement is 
reasonable, especially at higher biases, where the k = 1 
and 2 JQP-ph lines appear. Moreover, the simulation 
predicts the appearance of a feature at the intersection 
of the fc = 1 J-ph lines for Jl and J2 near Qo = —0.58 
and V « 374 /iV. A corresponding peak can also be 
seen in the data m Fig. ||(a). This feature is an ana- 
log of the 3e peakO but with Cooper pair tunneling re- 
placed by J-ph processes, and we refer to it as the 3e-ph 
peak. Our samples do not satisfy the condition i?c > |A 
for the appearance of the usual 3e cycle. The corre- 
sponding condition for the 3e-ph cycle in Sample 1 is 
i?c > |(A — hujs/^) = 115.6 /zeV which is just satisfied. 
A similar feature appears in the data (at Qo ~ —0.73 
and V « 440 /iV) near the intersection of a fc = 1 J-ph 
line for Jl and a fc = 2 J-ph line for J2. A corresponding 
peak appears at the same location in the simulation. 

Despite this level of agreement, some discrepancies ex- 
ist. The measured current remains substantial down to 
the 3e threshold line between 275-350 /iV, while the cal- 
culated current in Fig. |2|(b) nearly vanishes for all Qo 
below the 3e-ph peak. Furthermore, in the simulation 
the primary (k = 0) JQP peak is strongly suppressed 
at the intersections with J-ph transition lines; an exam- 
ple occurs in Fig. |(b) at Qo ~ -0.83 and V = 607 ^lY 
where the primary JQP peak for Jl intersects a fc = 2 line 
for J2. The standard JQP process would begin with the 
Cooper pair transition —1 ^ —3. Our simulations indi- 
cate that at this intersection the J-ph transition — 1 — > 1 
is allowed, and is followed rapidly by a 1 — > quasipar- 
ticle transition. The transition — s- — 1 is allowed but 
slow, so that the occupational probability of the n = — 1 
state is quite small (~ 0.01), and the JQP cycle is sup- 
pressed. There is, however, no clear sign of suppression 
of the primary JQP cycle in our data. 

We can account for these discrepancies to some ex- 
tent by considering PAT due to the environment, here 
assumed to be the 2DEG in the immediate vicinity of the 
SET. Since the effective environment temperature Xe can 
be much larger than the SET electron temperature, Q PAIi 
can be quite significant. Following Siewert and Sch6n,t2l 
we incorporate it into our calculation by adding the rate 
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FIG. 4: Simulated I vs. Qo for V = 607 /iV, with (dotted 
line) and without (dashed line) the efltects of PAT. Inset: 
measured / vs. Qo for Sample 1 at 1/ « 577 /xV (solid), and 
simulated 7 vs. Qo with (dotted) and without (dashed) PAT. 
Here Te = 2 K and i?e = 3 ^. 



to all quasiparticle tunneling events for which 5/ + 2 A > 
0. Fig. H(c) shows the results of a simulation including 
PAT with much improved agreement. Appreciable cur- 
rent appears at voltages around the 3e threshold lines. 

In addition, the suppression of the JQP peak at inter- 
sections with J-ph lines is largely lifted. This can be seen 
in greater detail in Fig. || which shows the calculated / 
vs. Qo for V = 607 ^V, both with and without PAT. 
The inclusion of PAT increases the primary JQP peaks 



for both Jl and J2. In the inset, we make a compari- 
son with data from Sample 1 at ^ « 577 ^V, a voltage 
below the cutoff for the JQP-ph cycle, so that only the 
JPQ peak is visible. The JQP peaks for both Jl and J2 
are near J-ph transition lines here, and in the absence 
of PAT the simulation underestimates the actual JQP 
current. Addition of PAT brings the simulation in much 
closer agreement with the data, although it does not cor- 
rectly predict the width of the JQP peak, a topic which 
will be addressed in a separate publication. 

In conclusion, we have investigated charge transport 
processes in an S-SET which is coupled to a 2DEG, 
so that the contact pad structures form microstrip res- 
onators. We have observed current-carrying cycles which 
involve tunneling of a Cooper pair and simultaneous 
photon emission, including JQP-ph and 3e-ph processes. 
Simulations performed using the orthodox theory agree 
well with the data, particularly when the effects of PAT 
due to the environment are included. In the simulations, 
suppression of the primary JQP in the vicinity of J-ph 
lines is partially lifted by the effects of PAT. While we 
do not observe the suppression, it may be detectable in 
samples for which the cental island of the SET is not so 
strongly coupled to an external environment. 
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